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Cold molecular gas in cooling flow clusters of galaxies 
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Abstract. Fhe results of a CO line survey in central cluster galaxies with cooling flows are presented. Cold 
molecular gas is detected with the IRAM 30m telescope, through CO(l-O) and CO(2-l) emission lines in 6-10 
among 32 galaxies. Fhe corresponding gas masses are between 3.10 8 and 4.10 10 Mq. Fhese results are in agreement 
with recent CO detections by |Edge (2001)| A strong correlation between the CO emission and the Ha luminosity 
is also confirmed. Cold gas exists in the center of cooling flow clusters and these detections may be interpreted as 
an evidence of the long searched very cold residual of the hot cooling gas. 



Key words. Galaxies, Clusters, Cooling flows, Molecular 
gas 

1. Introduction 



Studies of X-ray emission of hot intra-cluster medium 
(ICM) have pointed out the high density of this gas 
in the central regions of many clusters. The derived 
timescales for radiative cooling in the center is much 
smaller than the Hubble time, and the ICM is predicted 
to condense and flow towards the cluster center (see 
Fabian 1994 for a review). The X-ray spectra show 
evidence of cooler gas in the center, through central 
drops of temperature. But the fate of the cooled gas 
still remains uncertain. The duration of the cooling flows 
is thought to be a significant fraction of the cluster 
life-time, since cooling flows are quite frequent in clusters. 
Estimated cooling rates of the order of 100 M Q /yr and 
up to 1000 M /yr implied that enormous quantities of 
material should have accumulated (10 11 to 10 12 M Q in 
a fraction of a Hubble time). But no resulting cold gas 
has been detected in molecular form until recently. Many 
efforts have been expended to detect this gas in emis- 
sion or absorption, either in HI (Burns ct al. 1981 
Valentijn fe Giovanelli 1982| IShostak et al. 19831 

McNamara et al. 19901 IDwaraka nath et a l. 199511 o r 
in the CO molecule, see |Grabelsky fc Ulmer (1990 
McNamara et al. (1994)] |Antonucci et al . (199*3 



Uraine fc bupraz (Mi)} |0'Dea et al. Tfle 

intracluster medium is enriched in heavy elements with a 
metallicity of up to 0.3 solar making possible the forma- 
tion of CO molecules. The first detection of CO emission 



has been made in Perseus A by Lazareff et al. (1989) 
but the corresponding H2 is not strongly identified as 
coming from the cooling flow rather than from the galaxy 



itself. Recently, Edge (2001)| reported to have found CO 
line emission in the central galaxy of sixteen extreme 
cooling flow clusters. Starbursts that may appear as 
a consequence of the gas condensation must produce 
a lot of young and hot stars. But the observed stellar 
luminosities are not bright enough to account for the 
high mass deposition rates of cooling flows. Although 
Chandra and XMM-Newton observations lead to reduced 
rates, the cold molecular gas masses observed in some 
cluster cores remain a small fraction of the gas cooled 
along the flow. Th e ICM is probably multi-phase (e.g. 
Ferl and et al. 1994)l. A significant fracti on of gas might 
be so cold ( |Pfenniger fc Combes 1994] ) that it could 
correspond to the high concentration of dark matter 
in cl usters deduced from X-ray data and gr avitational 
arcs (IDurr et et al. 19941 IWu fc Ham mer 199 3|). Rec ently 
Lieu et al. (1996,1999)[ and |Mittaz et al (1998)] have 
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detected large quantities of gas at intermediate temper- 
ature of 5 • 10 5 K in 5 clusters with the EUVE satellite 
(Extreme Ultraviolet Explorer). Since this phase is quite 
transient, the mass flow implied would be much larger 
than that of the cooling flow itself. Other processes must 
be at work, such as heating by shocks, or mixing layer 
mechanisms at the interface between a cold and hot 
phase IjBonamente et al. 2001JI . Also the detection of the 
near-infrared quadrupolar emission line H 2 (1-0)S(1) in 
central cluster galaxies with cooling flows (and their non- 
detection in similar control galaxies without cooling flows, 
e.g. lFalcke et al. 1998(1 suppo rt the presence of mol ecular 
gas at tem perature of 2000K llJaffe fc Bremer 19971 Edge 
et al. 2002. IWilman et al. 2002^ 

In this paper, we present our search for CO lines in 32 
galaxies in the center of clusters, carried out in June and 
August 2001 with the IRAM 30m telescope. We have 
found 6 clear detections and 4 hints of CO lines. In the 
next section we describe the instrumental conditions of 
our observations and the data reduction. We then present 
results and cold Eras mass evaluations in section 3. In 
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Source 


Redshift 


RA (2000) 


Decl (2000) 


Central Freq. 


Exposure 


Central Freq. 


Exposure 














in C0(l-0) 


time (min) 


in CO(2-l) 


time (min) 


A85* 


0.05567 


00 


41 


50.4 


-09 18 11 


109.19 GHz 


152 






- 


Z235 


0.08300 


00 


43 


52.1 


+24 24 21 


106.43 GHz 


128 


212.87 


GHz 


128 


A262* 


0.01620 


01 


52 


46.5 


+36 09 07 


113.43 GHz 


160 






- 


A291* 


0.19590 


02 


01 


43.1 


-02 11 48 


96.38 GHz 


304 






- 


A496 


0.03281 


04 


33 


37.8 


-13 15 43 


111.60 GHz 


248 


223.21 


GHz 


256 


RXJ0439+05 


0.20800 


04 


39 


02.2 


+05 20 44 


95.42 GHz 


300 






- 


PKS0745-191 


0.10280 


07 


47 


31.3 


-19 17 40 


104.52 GHz 


304 


209.04 


GHz 


376 


A644* 


0.07040 


08 


17 


25.5 


-07 30 44 


107.69 GHz 


144 


215.37 


GHz 


88 


RXJ0821+07 


0.11000 


08 


21 


02.4 


+07 51 47 


103.84 GHz 


96 


207.69 


GHz 


96 


A646 


0.12680 


08 


22 


09.6 


+47 05 53 


102.30 GHz 


304 


204.59 


GHz 


312 


A780 


0.05384 


09 


18 


05.7 


-12 05 44 


109.38 GHz 


268 


218.76 


GHz 


352 


A 1068* 


0.13860 


10 


40 


44.5 


+39 57 11.1 


101.23 GHz 


180 






- 


A978 


0.05425 


10 


20 


26.5 


-06 31 36 


109.34 GHz 


96 


218.67 


GHz 


96 


A1668 


0.06368 


13 


03 


46.6 


+19 16 18 


108.37 GHz 


80 


216.73 


GHz 


80 


A 1795 


0.06326 


13 


48 


52.4 


+26 35 34 


108.41 GHz 


68 


216.82 


GHz 


129 


A2029 


0.07795 


15 


10 


56.1 


+05 44 41 


106.93 GHz 


112 


213.86 


GHz 


168 


MKW3s 


0.04531 


15 


21 


51.9 


+07 42 32 


110.27 GHz 


144 


220.54 


GHz 


144 


A2146 


0.23370 


15 


56 


13.8 


+66 20 55 


93.43 GHz 


272 






- 


A2142* 


0.09037 


15 


58 


20.0 


+27 14 02 


105.71 GHz 


128 






- 


A2147 


0.03532 


16 


02 


17.0 


+15 58 28 


111.33 GHz 


152 


222.67 


GHz 


152 


A2151 


0.03533 


16 


04 


35.8 


+17 43 18 


111.33 GHz 


152 


222.67 


GHz 


152 


A2199 


0.03035 


16 


28 


38.5 


+39 33 06 


111.87 GHz 


128 


223.74 


GHz 


120 


Z8193* 


0.18290 


17 


17 


19.2 


+42 27 00 


97.44 GHz 


304 








A2261* 


0.2240 


17 


22 


27.1 


+32 07 58 


94.17 GHz 


256 








A2319* 


0.05459 


19 


21 


10.0 


+43 56 44 


109.30 GHz 


136 


218.60 


GHz 


88 


CygA 


0.05607 


19 


59 


28.3 


+40 44 02 


109.15 GHz 


228 


218.29 


GHz 


232 


A2462 


0.07437 


22 


39 


11.4 


-17 20 28 


107.29 GHz 


144 


214.58 


GHz 


144 


A2597 


0.08520 


23 


25 


19.8 


-12 07 26 


106.22 GHz 


144 


212.43 


GHz 


144 


A2626 


0.05490 


23 


26 


30.6 


+21 08 50 


109.27 GHz 


136 


218.54 


GHz 


140 


A2634 


0.03022 


23 


38 


29.5 


+27 01 56 


111.89 GHz 


144 


223.77 


GHz 


144 


A2657 


0.04023 


23 


44 


57.4 


+09 11 34 


110.81 GHz 


208 


221.62 


GHz 


208 


A2665* 


0.05610 


23 


50 


50.6 


+06 09 00 


109.14 GHz 


192 


218.29 


GHz 


192 



Table 1. This table presents the sample of cooling flow clusters of galaxies observed with the IRAM 30m telescope. 
Central frequencies of the CO(1-0) and CO(2-l) lines observed, as the exposure time are indicated for each galaxy. 
Sources observed in the second run (August) are with indicated by a star. The others were observed during the first 
run (July). Several sources were not observed in CO(2-l) since the redshifted J=2-l transition lines were out of the 
30m telescope receiver's band. 



section 4 and 5 we discuss the possible significations 
of such large gas quantities when they are present and 
compare these measurements with other wavelength 
observations. 



2. Observations and data reduction 

The sample of sources was selected according to several 
criteria. First, we wanted to observe galaxies with im- 
portant cooling flows, so we chose high deposition rates 
galaxies with M around or greater than 100 M©/yr, 
see |Peres et al. (1998)| |White et al. (1997)| though these 
rates are certainly overestimated. Three non-cooling flow 
clusters (Abell 1668, Abell 1704 and Abell 2256) have also 
been observed and not detected in CO with the 30m tele- 
scope. It is possible that the large gas flow produces mas- 
sive stars ionizing the gas. The gas might also be cooling 
in ionizing shocks (optically luminous). Thus, the pres- 



nied by Ha emission as suggested in Edge (2001) Sources 



were then selected according to their H Q luminosity when 
available (high luminosity of about 10 42 erg.s" 1 from 



Crawford et al. (1999)| |Owen et al. (1995)| ). The sample 
contains only relatively low-redshift cD galaxies (z < 
0.25), for the sake of sensitivity. We gather data at other 
wavelengths, such as the far infra-red, when available, to 
be able to compare gas and dust emission. All observing 
parameters are summarized in Tabled Observations were 
achieved with the IRAM 30m millimeter-wave telescope 
at Pico Veleta, Spain in June and August 2001 in good 
weather conditions. We used four receivers simultaneously, 
centered two on the CO(1-0) and two on the CO(2-l) lines 
at 115 GHz and 230 GHz. The beam of the telescope at 
these two frequencies is 22" and 13" respectively. Two 
backends were provided by the autocorrelator, with a 1.25 
MHz resolution on a 600 MHz band width. The two other 
backends were the two 512MHz wide 1MHz filter-banks. 
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Source 


Line 


Peak 


Rms 


Line 


Line 


Line 


Ico 






mK 


mK 


detection 


position (km/s) 


width (km/s) 


K.km/s 



A85* 


CO 


1-0) 




0.9 


no 




300* 


<0.35 


Z235 


CO 


1-0) 


1.9±0.2 


0.8 


hint 


-258±48 


318±103 


0.66±0.19 




CO 


2-1) 




1.7 


no 




300* 


<0.65 


A262* 


CO 


1-0) 


2.9±0.5 


0.8 


yes 


31±24 


346±47 


1.05±0.14 


A291* 


CO 


1-0) 




0.9 


no 




300* 


<0.35 


A496 


CO 


;i-o) 


1.5±0.5 


0.7 


hint 


382±39 


311±73 


0.49±0.12 




CO 


[2-1) 


3.0±0.7 


1.1 


hint 


114±29 


249±52 


0.80±0.18 


RXJ0439+05 


CO 


1-0) 




0.72 


no 




300* 


<0.28 


PKS0745-191 


CO 


1-0) 


2.0±0.3 


0.6 


yes 


18±29 


221±58 


0.47±0.11 




CO 


2-1) 


11.2±0.8 


2.1 


yes 


-45±15 


215±41 


2.57±0.38 


A644* 


CO 


1-0) 


1.1±0.4 


0.7 


no 


-10±55 


260±77 


<0.27 




CO 


2-1) 




3.1 


no 




300* 


<1.12 


RXJ0821+07 


CO 


1-0) 


8.9±0.6 


1.1 


yes 


-2± 8 


135±21 


1.28±0.16 




CO 


2-1) 


9.9±1.5 


2.9 


yes 


5±34 


270±81 


2.86±0.60 


A646 


CO 


1-0) 


1.5±0.1 


0.5 


yes 


105±39 


376±121 


0.62±0.15 




CO 


2-1) 


1.9±1.0 


1.6 


no 


46±54 


346±123 


<0.62 


A780 


CO 


1-0) 


1.8±0.5 


0.6 


yes 


219±31 


439±68 


0.80±0.11 




CO 


2-1) 




1.8 


no 




300* 


<0.69 


A978 


CO 


1-0) 


1.8±0.7 


1.0 


no 


-128±49 


221±87 


<0.39 




CO 


2-1) 




2.1 


no 




300* 


<0.80 


A1068* 


CO 


1-0) 


10.1±0.3 


0.5 


yes 


-45±4 


249±10 


2.66±0.1 


A1668 


CO 


1-0) 




1.2 


no 




300* 


<0.46 




CO 


2-1) 


2.4±1.6 


1.7 


no 


9±64 


281±125 


<0.65 


A1795 


CO 


1-0) 


3.4±0.7 


0.9 


yes 


-190±28 


405±56 


1.47±0.19 




CO 


2-1) 


6.1±1.0 


1.8 


yes 


-128±31 


500±65 


3.26±0.39 



Table 2. Summary of observational data for the two runs. Lines characteristics are presented, spectra are shown on 
Fig-El Ico (column 8) were evaluated from a gaussian fit of the CO(1-0) and CO(2-l) lines for detections and hints of 
detections. Ico upper limits were evaluated by equation [3 for non-detections. A detection is asserted in one transition 
line when the peak of the gaussian fit is above three times the rms (in 55 km/s channels) and a hint of detection 
when it is between two and three times the rms. To claim a cold molecular gas detection, we require a detection in 
both transition lines. A possible detection was claimed when was present a detection in one transition line or a hint 
of detection in both lines, unless the line appear clearly in one transition only (Abell 262, Abell 1068, Zw8193). 



~ 650km/s at 1.3mm. In addition, we used the 4 MHz 
resolution filter-bank, providing a 1 GHz band width, im- 
portant for the 1.3mm receivers (since it corresponds to 
1300km/s bandwidth also). Given the uncertainty in the 
central velocity of the CO line (some optically measured 
velocities being systematically displaced with respect to 
the galaxy systemic velocity), the expected width of a cD 
galaxy, and the required baseline to eliminate sinusoidal 
fluctuations, this wide band is necessary. The signals are 
expressed in main beam temperatures, since the sources 
are not expected to be extended and homogeneous. The 
main-beam efficiency of the 30m is : 



Vmb = Tl/T mh and r] mh = B eff /F eff 



(1) 



with ?7 mb = 0.75/0.95 at 115 GHz and 0.52/0.91 at 230 GHz 
(cf IRAM-30m site http:/ /www.iram.es/). The data were 
reduced with the CLASS package; spiky channels and bad 
scans were excised. After averaging all the raw spectra for 
each line of each source, linear baselines were subtracted 
and the spectra were Hanning smoothed. Assuming that 
a good sampling of the line requires at least five points, 
and assuming a typical 300 km/s line width, the data 



and 2.6mm, to gain more than a factor 3 in signal to 
noise. The CO emission lines were fitted with gaussian 
profiles through the CLASS package. For non detections, 
CO intensity Ico upper limits at 3<r were evaluated as in 
McNamara et al. (1994)] by : 



I 



CO 



3 x rms x W Une x ( 



W s , 



a/2 



(2) 



where rms is the noise level, computed in the channel 
width W sm (55 km/s) and Wu ne is the expected line band- 
width, typically 300 km/s. The typical integration time 
was two hours for all galaxies. CO intensities were ob- 
tained by averaging the data of at least two different days, 
to prevent systematic effects (baseline ripples...). The de- 
tection criterion is the CO emission line maximum is at 
least three times the rms, in 55km/s channels. When the 
signal to noise ratio was between 2 and 3 rms we con- 
cluded to a hint of detection. Below 2 rms, upper lim- 
its were computed for the CO intensity. Results are pre- 
sented in Table and Table 01 After this first analysis at 
each wavelength, we compare Ico(i-o) an( A Ico(2-i)- A 
cold gas detection is claimed when CO is detected in the 
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Source 


Line 


Peak 


Rms 


Line 


Line 


Line 


Ico 






mK 


mK 


detection 


position (km/s) 


width (km/s) 


K.km/s 


A2029 


CO(l-O) 


1.1±0.5 


0.9 


no 


-68±87 


445±168 


<0.35 




CO(2-l) 




1.6 


no 




300* 


<0.62 


MKW3s 


CO(l-O) 


1.3±0.4 


0.8 


no 


-133±60 


347±100 


<0.31 




CO(2-l) 




1.0 


no 




300* 


<0.39 


A2146 


CO(l-O) 




0.5 


no 


_ 


300* 


<0.19 


A2142* 


CO(l-O) 




0.7 


no 




300* 


<0.27 


A2147 


CO(l-O) 




1.1 


no 




300* 


<0.42 




CO(2-l) 




1.4 


no 




300* 


<0.54 


A2151 


CO(l-O) 


2.0±0.3 


0.8 


hint 


341±47 


133±95 


0.29±0.14 




CO(2-l) 




1.5 


no 




300* 


<0.58 


A2199 


CO(l-O) 


1.6±0.4 


0.9 


no 


109±55 


395±123 


<0.35 




CO(2-l) 


3.1±0.6 


1.2 


hint 


265±28 


204±61 


0.67±0.18 


Z8193* 


CO(l-O) 


2.1±0.4 


0.4 


yes 


-14±20 


242±65 


0.55±0.1 


A2261* 


CO(l-O) 




0.8 


no 




300* 


<0.31 


A2319* 


CO(l-O) 




0.7 


no 




300* 


<0.27 




CO(2-l) 




3.4 


no 




300* 


<1.3 


CygA 


CO (1-0) 


3.6±0.3 


1.4 


hint 


203±24 


153±46 


0.59±0.18 




CO(2-l) 


2.1±0.5 


1.4 


no 


-83±65 


336±207 


<0.54 


A2462 


CO(l-O) 




0.9 


no 




300* 


<0.35 




CO(2-l) 




2.7 


no 




300* 


<1.0 


A2597 


CO(l-O) 


1.1±0.3 


0.8 


no 


-10±55 


260±77 


0.29±0.11 




CO(2-l) 




2.3 


no 




300* 


<0.89 


A2626 


CO(l-O) 


1.7±0.2 


0.8 


hint 


-161±42 


293±103 


0.52±0.15 




CO(2-l) 




1.4 


no 




300* 


<0.54 


A2634 


CO(l-O) 


2.9±0.4 


1.1 


hint 


488±62 


459±137 


1.44±0.38 




CO(2-l) 




2.2 


no 




300* 


<0.85 


A2657 


CO(l-O) 


3.0±0.03 


0.8 


yes 


147±13 


69±50 


0.22±0.080 




CO(2-l) 




0.8 


no 




300* 


<0.31 


A2665* 


CO(l-O) 


1.5±0.2 


0.8 


no 


-273±40 


193±82 


<0.31 




CO 2-1 




1.6 


no 




300* 


<0.62 



Table 3. Summary of observational data for the two runs. Lines characteristics (Table El continuation) . 



detections are defined by a detection in one line or by a 
hint of detection in both lines. Otherwise, we concluded to 
a non detection. For Abell 262, Abell 1068 and Zw8193, 
without CO(2-l) data, a detection is claimed because of 
a clear line detection in CO(1-0). With these criteria, we 
claim 6 detections, 4 possible detections, and 22 non de- 



tections. Detections made by Edge (2001) for Abell 1068 
and RXJ0821+07 are confirmed, with a good agreement 
of the derived hydrogen molecular gas masses. CO (1-0) 
emission line is also confirmed in Abell 262, but molecular 
gas mass deduced here is twice lower than in Edge (2001) 



width. Three new values of M gas are found : Abell 646 for 



which cold gas mass is in agreement with E dge (2001)| up- 
per limit, Abell 1795 for which a large line width is found 
and consequently a molecular hydrogen mass higher than 



the upper limit deduced in|Edge (2001]] and PKS0745-19 



that was observed at a wrong frequency by Edge (2001) 



CO intensities always correspond to areas deduced 
from gaussian fits for detections and possible detections 
when a gaussian fit was possible, but for the 22 no- 
detections, the I co are evaluated with Formula (2). Table 
0] show gas mass estimates as well as X-ray, optical, IR 



3. Results 

3.1. H2 mass evaluation from CO observations 

Since cold H 2 is a symmetric molecule, the best tracer of 
cold molecular gas is the CO lines, from the most abun- 
dant molecule after H 2 : CO/H 2 ~ 6.10~ 5 . From standard 
(and empirical) calibrations, it is possible to deduce the 
interstellar H 2 content from the integrated CO intensity 
I co (K.km/s) : 



That comes from the fact we identify a line with a smaller ^ co 



T mh (CO)dV 



(3) 



where T m ^(CO) is the main beam antenna temperature, 
obtained for the CO(1-0) line (cf Sect. 2). Although the 
typical molecular cloud is optically thick in the first CO 
lines of the J ladder, the proportionality factor between 
the column density of molecular gas and integrated inten- 
sity is justified, since the observed signal is the emission 
sum of many clouds in the beam, and these clouds have 
a small filling factor, when spatial and velocity volume 
is considered. We adopt here the conversion factor com- 
monly used for N(H 2 ) in molecule/cm 2 unit : 
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From this equation the mass of molecular hydrogen, con- 
tained in one beam, in M Q is: 



(1999) 



M(H 2 ) = 2.95 lCr 19 I co 6 2 D 2 N ^ 



I 



CO 



(5) 



where Ico is the integrated intensity in K.km/s, 8 is the 
beamsize of the telescope in arcseconds, and D is the dis- 
tance of the galaxy in Mpc, determined with a Hubble 
constant Hq—70 km/s/Mpc. This converting factor is in 
agreement with the one used in Edge (2001) and in pre- 
vious CO observations of cluster cores but gas mass es- 
timations are slightly lower here because we took Hq=70 
km/s/Mpc. It is important to notice the Ico/N(H.2) con- 
version factor has been first calibrated in the solar neigh- 
borhood. To use this value implies intra-cluster medium 
is assumed to behave like the Galactic interstellar mat- 
ter near the sun, in particular with the solar metallicity. 
But the intracluster medium has subsolar metallicity. So 
this conversion factor is likely to underestimate the mass 
of molecular hydrogen. A standard factor 1.36 taking into 
account He contribution is also used in the gas mass esti- 



optical lines, see Crawford et al 
PKS 0745-191 i s a 0.1028 re dshifted galaxy already 
observed in CO by Edge (2001)| In present observations, 
a CO(1-0) line is detected three times above the rms. A 
simultaneous detection in the CO(2-l) band with a signal 
to noise better then five seems to confirm the presence 
of molecular gas. This galaxy is supposed to contain a 
large cooling flow with mass deposition rates around 1000 
M /yr according to |Peres et al. (1998)| |Allen (2000)] 
Strong optical emission lines have also been detected in 
PKS 0745-191. This galaxy is the site of an important ex- 
citation mechanism. Besides it is a powerful radio source 
with an amorphous and filamentary morphology. 

Recently, 



sec 1 



Baum fc O'Dea (1991) 



Donahue et al. (2000) 



mation : M ffQS =1.36M(H 2 ). The derived M gas values, dis- CO by |Edge (2001) 
played in Table H are between ~ 10 s and ~ 10 10 M . 
These masses, together with upper limits, are plotted on 
Fig. H for all sources, as a function of redshift. The curve 
represents the M sas detection limit of the IRAM 30m tele- 
scope for the CO(1-0) line, assuming a typical 300 km/s 
lincwith and a temperature detection limit of 0.5 mK, in 
agreement with our noise level. 



have mapped kpc-size filaments in vibrationally-excited 
H2 in the cores of galaxies centers of cooling flows, like 
PKS 0745-191, with high spatial resolution. They have 
also found dust lanes which are optically thick to 1.6/im 
emission. These dust lanes are confined to the central few 
kpes. 

The cD galaxy RX J0821+07 has been detected in 
This relatively easy detection is 



confirmed here in both wavelength and CO intensities 
are compatible with previous ones. Optical images 
taken with the AAT and Hubble Space Telescope by 
Bayer-Kim et al. (2002)| show that the central galaxy 



is embedded in a luminous and extended line-emitting 
nebula coincident with a bright excess of X-ray emission 
imaged by Chandra. 



Abell 1068 detected by Edge (2001) is confirmed here. 
The H2 mass deduced is the highest we found among 
detections. This galaxy show strong optical lines and 
large dust mass in comparison with other detections. It 
is also a powerful IRAS source with a 650 mJy flux at 60 
/mi. 

Abell 1795 has been observed by 



Braine & Dupraz (1994) who did not detect CO line. 



Fig. 1. M gas in 1O 8 M deduced from CO observations, 
versus the galaxy redshift. Filled stars are detections, filled 
square hints of detection and arrows are upper limits. The 
dashed line represents the molecular mass limit than can 
be deduced from CO observations with the IRAM 30m 
telescope, in 2h integration time. 



3.2. CO detections 



Edge (2001) found a marginal detection. This detection is 
confirmed here in the two bands. But line widths deduced 
from gaussian fits are quite different in the two wave- 
lengths, so we cannot exclude the possibility of very high 
velocity molecular clouds. This galaxy is known as a radio 

An optical filament has 



David et al. (1993) 



Abell 262 has also been detected by Edge (2001) 



Current values are compatible with these measurements 
and confirm this detection even if CO (2-1) data were 
not good enough in the second run of August 2001 to 
detect a line. This galaxy contains a central radio source 



source see 

been detected in H a by Cowie et al. (1983) According to 
them, some of the filaments observed in Abell 1795 seem 
to be concentrated and coming from the galaxy whereas 
fainter extended filaments are surrounding the galaxy. 
The question about their origin and their link with the 
cooling flows is not clearly determined. Mapping cold gas 
will allow to better understand the spatial structure of 
the cooling material and to know if the CO is along the 
filaments or in the galaxy, Salome & Combes (2003, in 
prep). 

Zwicky 8193 is a strong optical line emitter. Only CO(l- 
0) was observed here. The molecular gas mass deduced 
agree with the value derived in |Edge (2001)1 Zwicky 8193 
is a complex system and we refer to lEdge (2001)1 and 
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3.3. Hints of detection 



4.2. Optically thick cold clouds 



We consider 4 galaxies of our sample to be possible CO 
emitters according to criteria denned above. Nevertheless 
CO emission lines here are fainter than the previous ones, 
with values reaching half a K.km/s. 

Abell 496 is a cD galaxy. A possible line is seen in 
the two bands, but signal to noise ratio between 2 and 
3 is not sufficient to claim a detection. Much time has 
been dedicated to this radio source, see for example 
Peres et al. (1998) to deduce a small upper limit of H2 
mass. Faint optical H Q line have been observed in this 



galaxy also emitting in X-ray, see David et al. (1993) 



Nevertheless we deduce here a new upper limit in molec- 
ular gas mass. 

We also assert a possible detection in Abell 646, even 
if no CO(2-l) line is seen, because of the clear shape 
of the CO(l-O) line detected just above three time the 
rms. Moreover, Edge (2001) asserted to have a marginal 



detection of this galaxy. So it would be interesting to 
confirm this detection. 

Abell 780 (Hydra A) is a very powerful radio source 
that had already been observed through millimetric 
wavelength, see for example O'Dea et al. (1994b) This 
much studied source is here at the limit of detection 
in CO(l-O) and not seen in CO(2-l). M gas upper limit 
deduced is in agreement with the evaluation made by 
Edge (2001)[ but no clear detection can be claimed. 
The cD Abell 2657 galaxy with optical emission lines 
was not detected in CO(2-l). Faint possible CO(l-O) line 
is present and a new upper limit in M gas is derived, but 
more observations are required. 



3.4. Upper limits 

A large number of the selected galaxies are not detected. It 
is possible, these cooling flow clusters of galaxies contains 
cold gas with lines still too weak to be detectable with the 
actual IRAM 30m telescope sensitivity. 

4. Discussion 

4.1. Origin of the cold gas 

Recent X-ray observations by Chandra, see 
IFabian 20021 |Voigt et al. (2002)| and X MM-Newton 
UPeterson et al. 200*TI ITamura et al. 20 01) have con- 
firmed the presence of radial gradients in temperature in 
the cores of several clusters of galaxies. Even if results 
from the high spectral resolution Reflection Grating 
Spectrometer (RGS) on XMM-Newton do not show 
evidence (from Fe XVII) for gas cooling at temperature 
lower than 1-2 kev, millimetric emission of a cold gas 
component is detected in the center of several galaxy 
clusters. Added to recent data from Edge (2001)[ new 
detections of molecular gas in cooling flow galaxies is of 
great interest. But questions persist on the origin of this 



In an optically thick medium which is the case here, the 
CO(2-1)/CO(1-0) ratio should be about or less than one, 
if we assume the same excitation temperature for the 
two CO energy levels. On Fig. is plotted the CO (2- 
1) versus CO(1-0) intensity (in K.km.s -1 ) for the galax- 
ies observed here. The straight line indicates their equal- 
ity. CO(2-l) intensities have been multiplied by a beam 
correcting factor ~4 and by the relative beam efficien- 
cies (0.52/0. 91)/(0.75/0.95)=0.72 to be compared to the 
CO(1-0) ones. The preliminary plot indicates that the 
CO(2-l) line is in fact lower than the CO(1-0) one; this is 
in general the case for sub-thermally excited gas, in nearby 
galaxies fe.g. lBraine & Combes 1992|l . The CO lines ratio 
are consistent with an optically thick gas. The medium 
considered here is certainly far more complex, probably 
inhomogeneous and multi-phase. It might be a mixing of 
diffuse gas and denser clumps, and the diffuse medium 
might be dominating the emission, while thick and small 
clouds could enclose a larger quantity of hydrogen mass 
than estimated. 




CO(l-O) (K.km/s) 



Fig. 2. CO(2-l) versus CO(1-0) corrected from the effect 
of different beam sizes. The straight line corresponds to 
line emission equality. 



4.3. Origin of cold gas 

If the molecular gas is formed by the cooling flow, 
we should see a correlation between the detected 
cold gas masses and the X-ray determined mass de- 

But 



position rate, as shown by Edge (2001) But there 
is a quite large dispersion in the Mx values be- 
cause of the different methods used in the literature, 
see |Grabelsky fc Ulmer (1990)| |Bregman et al. (1990 



White et al. (1997)1 l Peres et al - (1998)1 or |Allen (2000J 



To test this, we have compared cold gas masses founc 
here with mass deposition rates evaluated thanks to an 
Einstein Observatory X-ray image deprojection analysis 
made by |White et al. (1997)| see Fig. EH The comparison 
is also done with M issue from a ROSAT observatory 



spatial analysis by Peres et al. (1998) as shown on Fig. 
0] These two samples of mass deposition rates, evaluated 
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Source 


M{H 2 ) 


Mdust 


L(H a ) 


M« 


(i) 

cool 


M< 2 > 


(2) 
cool 


Flux(1.4GHz) 




(M Q ) 


(M Q ) 


(erg/s) 


(M /yr) 


(kpc) 


(M Q /yr) 


(kpc) 


(mJy) 


A262 


2.3±0.3xlO s 


8.7xlO b 


6.0xl0 39 


27 


104 


10 


67 


131 


pks0745-191 


4.0±0.9xl0 9 


3.5xl0 7 


1.4xl0 42 


1038 


214 


579 


177 


2370 


rxj0821+07 


1.3±0.2xl0 10 


4.2xl0 7 


3.0xl0 41 


- 


- 


- 


- 


- 


A1068 


4.2±0.2xl0 10 


1.4x10 s 


1.7xl0 42 


- 


- 


- 


- 


- 


A 1795 


4.8±0.6xl0 9 


6.7xl0 7 


l.lxlO 41 


381 


177 


321 


181 


930 


Z8193 


1.5±0.3xl0 10 


3.8xl0 7 


1.5xl0 42 


- 


- 


- 


- 


- 


A496 


4.3±1.0xlO s 


- 


3.4xl0 42 


95 


110 


134 


138 


- 


A646 


7.9±0.2xl0 9 


2.5xl0 7 


1.6xl0 41 


- 


- 


- 


- 


- 


A780 


2.0±0.3xl0 9 


3.0xl0 7 


1.6xl0 41 


262 


162 


222 


170 


40800 


A2657 


2.9±0.1xl0 8 


- 


- 


- 


- 


44 


101 


- 


A85 


< 8.8 x10 s 


- 


- 


107 


93 


108 


131 


58 


Z235 


< 2.5 xlO 9 


- 


- 


- 


- 


- 


- 


- 


A291 


< l.lxlO 10 


84.9 xlO 7 


4.6xl0 41 


- 


- 


- 


- 


- 


A644 


< l.lxlO 9 


- 


- 


189 


141 


136 


Ill 


- 


rxj0439+05 


< 9.8xl0 9 


5.1xl0 7 


l.lxlO 42 


- 


- 


- 


- 


- 


A978 


< 9.3xl0 8 


- 


2.0xl0 40 


- 


- 


- 


- 


- 


A1668 


< 1.5xl0 9 


- 


1.2 xlO 41 


- 


- 


- 


- 


- 


A2029 


< 1.7xl0 9 


- 


8.0 xlO 39 


555 


186 


431 


192 


550 


MKW3s 


< 5.2x10 s 


- 


5.0xl0 39 


175 


171 


132 


158 


- 


A2146 


< 8.6xl0 9 


8.8xl0 7 


1.4xl0 42 


- 


- 


- 


- 


- 


A2142 


< 1.8xl0 9 


- 


- 


350 


150 


369 


172 


- 


A2147 


< 4.3 xlO 8 


_ 


7.0xl0 39 


_ 


_ 


- 


_ 


_ 


A2151 


< 3.1x10 s 




5.8xl0 40 






166 


146 




A2199 


< 2.6 x10 s 




3.5xl0 40 


154 


143 


97 


124 


3700 


A2261 


< 1.3xl0 10 








20 


53 






A2319 


< 6.6x10 s 




l.OxlO 41 












CygA 


< 1.4xl0 9 




6.5xl0 42 


244 


135 


242 


167 




A2462 


< 1.6xl0 9 




5.8xl0 40 












A2597 


< 1.8xl0 9 


8.4xl0 7 


5.2xl0 41 


271 


152 






1880 


A2626 


< 7.6 x10 s 




3.3xl0 40 






53 


114 




A2634 


< 3.2x10 s 




3.7xl0 40 










7657 


A2665 


< 7.9x10 s 




6.0xl0 39 













Table 4. Derived parameters of the observed sources. Optical line luminosities are from |Crawford et al. (1999) Mass 
deposition rates and cooling radius are from (1) Peres et al. (1998) and (2) White et al. (1997). Dust masses are 
evaluated from 60/im data compiled in Edge (2001), assuming T £ i„ st =35K. 



the largest number of sources in common with the clus- 
ters observed here in CO. For a correlation trend to be 
relevant, the aim was (i) to have a high number of sources 
observed in both CO and X-ray, with regards to the faint 
detection level in CO and (ii) to compare M(H2) with 
M derived from one method only (with the same criteria 
for all sources). We can see a trend of correlation do ap- 
pear, even if there are very few data points. This confirms 
the relation between the mass of the cold component and 
the mass deposition rate already noticed in Edge (2001). 
Galaxies for which measurements have been possible lie 
close to M gas = l%xMxlGyr (large symbols). Then, as- 
suming that simple models of a multiphase flow would 
lead to an integrated mass deposition profile of the form 
M(<r)ocr Q , with a ~1. We have re-evaluated what would 
be the M inside the 30m telescope radius with a simple 
scaling by the cooling radius to the CO radius ratio. The 
correlation still appear but the cold gas masses detected 
are now close to M gos =Mx 10%x lGyr (small symbols 
with gray background) . These mass deposition rates have 



probably been overestimated by about a factor 5-10, see 
|McNamara et al. (2002)| as suggested by the recent X-ray 
observations by Chandra and XMM-Newton (e.g. Abell 
1795 in lFabian 20021 Abell 2199 in| Johnstone et al. (2002) 
or Abell 496 in |Dupke & White (2003)| ). Taking into ac- 
count the uncertainty on the conversion factor between 
H2 and CO, as discussed above, the correlation is in ac- 
cordance with a cooling scenario in which hot gas lead 
to cold substructures at rate deduced by X-ray observa- 
tions and detected here in CO (for an assumed age of the 
cooling is a few Gyr in the central regions). 



But many galaxies in cooling flow clusters observed 
here, do not show CO emission lines. Given the faint emis- 
sion temperature, it is possible that the cold gas is present 
but its radiation is below the detection limit. An alter- 
native is the gas is not cooled identically in all clusters 
of galaxies centers, depending on the environment of the 
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cold molecular gas mass detected, and for stronger ra- 
dio sources, M gas decreases when the radio power in- 
creases. This is consistent with a self-regulated heat- 
ing model powered by a central AGN (as suggested by 
Boh ringer et al. (2001)[ ) . But heating the ICM is certainly 



due to the radio lobes expansion whose energy is not only 
linked to the radio power. 



dM x /dt (M e /yr) 

Fig. 3. M gas from CO observations with respect to M the 
mass deposition rates deduced from Einstein X-ray data 



by |White et al. (19977] Straight lines from right to left 



are for M 



gas- 



ixMxlGyr, M gas = 10%xMx lGyr and 



=MxlGyr 




dKL/dt (M s /yr) 



Fig. 4. M, 

mass deposition rates deduced from ROSAT X-ray data 



g as from CO observations with respect to M the 



by |Peres et al. (1998)| 



4.4. Heating by AGN 

Many studies developed recently are taking into account 
heating mechanisms in cooling flows that could slow 
down the cooling and eventually stop it, which might 
explain the lack of CO emission in some of the galax- 
ies observed here. Important absorption could also hide 
the gas lying in a colder X-ray phase. Physical con- 
ditions in the central regions are certainly very com- 
plex, and simple cooling appears to be insufficient to 
explain multi- wavelength observations. Chandra images 
from Fabian et al. (2002)| or Johnstone et al. (2002) have 
pointed out holes coincident with radio lobes and cold 
fronts showing the interaction between the radio source 
and the intra-cluster medium (e.g. in Hydra A, the Perseus 
cluster, Abell 1795, Abell 2199, the Virgo cluster). It 
seems that the radio source and jets could heat the gas 
with shocks and significantly decrease the cooling rates 
IjDavid et al~200Tl |Briiggen et al. (2002)] ). Besides, 71 % 
of central cD galaxy in cooling flow clusters show a strong 
radio activity compared to 23% for non-cooling flow clus- 
ter cDs (Ball et al. 1993|l . No correlation has been found 
here between molecular gas masses and radio power at 
1.4 GHz (see fig. El- Nevertheless, it seems that for faint 



p ..4r,H, ( ff Hz '> 

Fig. 5. Mqas versus radio power at 1.4 GHz. There is no 
clear correlation between molecular gas masses and radio 
flux. 



4.5. Correlations with H a and Md ust 

The molecular hydrogen mass appears to be correlated 
with the amount of ionized gas (see Fig. EJ. Excitation 
processes leading to the gas emission could be the same for 
the two phases. The Ha emission could come from shocked 
cooling gas. The gas may also have been ionized from mas- 
sive stars born in a starburst triggered by the cooling flow. 
The best linear fit between the two components is plot- 



ted including previous detections of Edge (2001) (Figure 
El . Heating by a young star population is often suggested 
(see |Johnstone et al. (2002)| | . In that sense, gradients of 
metallicity deduced from Chandra observations could be 
explained by SN la injection of metal in the central galaxy 
(with the condition of some exchange of the gas at differ- 
ent radius, and so a possible mixing of different phases of 
the gas if they are present). The cold gas detected here 
might be a reservoir available for such a star formation 
process. 

The intracluster gas should be depleted in dust, at a 
given metallicity, since in the ICM environment and its 
physical conditions, the sputtering time of dust is much 
shorter than the dynamical time. The gas coming from a 
cooling flow, already at low metallicity, is thought to have 
a large relative depletion in dust. Therefore, the expected 
ratio between the CO measured gas content and the dust 
content from its submillimeter or far-infrared emission is 
large. Dust masses, derived from IRAS are evaluated for 
two assumed dust temperatures Tdust = 35 and 40 K. 
Theses masses are compared to cold gas masses on Fig.0 
It is important to notice how much dust mass highly de- 
pends on dust temperature. The gas-to-dust ratio for both 
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^*11Z2 4AI668 J, 

A2665 <f "3>^ A-4zB l.,r, c 

1 ^^"^ JA496A2J19 I 

A A2634AZJ51 J, 



L(Ha) (lO^erg.s- 1 ) 

Fig. 6. Mgas deduced from our CO observations versus 
H a . The straight line here is a fit to the detection and 
hint of detection points. Data from |Edge (2001)1 are plot- 
ted. Diamond shaped are the detections and triangles are 
the upper limits. Upper limits were not taken into account 
in the fit. There is a clear correlation between H Q lumi- 
nosity and cold gas masses, in agreement with previous 
detections. 









|J646^f m04 j|p«6 


#1795 





only IRAS data have been used here, tracing the warm 
dust. Significant amount of cold dust might be present 
as suggested by JCMT SCUBA detections in Abell 1835, 
Abell 2390 (Edge et al. 1999). More longer wavelength ob- 
servations at 850 (im tracing this cool dust would be of 
great interest. Nevertheless the mass to dust ratio found 
here are not incompatible with a cooling flow origin of the 
molecular gas. Besides, there is a trend of correlation be- 
tween cold gas masses and dust masses, but with a large 
dispersion. Infrared emission might be tied to star forma- 
tion. In that sense, Fig. [7| could also be interpreted as a 
possible correlation between gas content and star forma- 
tion, probably very active as we have seen previously. 

4.6. Excited H 2 detected in NIR 

Some of our sample galaxies have strong molecu- 
lar hydrogen emission in the 2 /im 1-0 S(l) line, 
HElston fc Maloney 1992| l994). This excited gas is 
thought to be associated to the cooling flow, since it is 
not detected in non-cooling flow galaxies of similar-type 
IjJaffe fc Bremer 1997|l . They reveal dusty nebular fila- 
ments, very similar to those detected in early type galaxies 
in small groups e.g. |Goudfrooij fc T rinchicr i~1998[ and in 
interacting gas-rich galaxies. The filaments are extended 
over kiloparsecs, and their heating source is not known. 
Recently, Donahue et al. (2000) have mapped the kpc- 



size filaments in vibrationally-excited H2 in the cores of 
galaxies centers of Abell 2597 and PKS 0745-191 with 
high spatial resolution. They have also found dust lanes 
which are optically thick to 1.6/im emission, confined to 
the central few kpes. Excited H2 produced directly by 
the cooling flow seems difficult, since H2 is much too 
luminous, by at least 2 orders of magnitude. It cannot be 
AGN photoionization or fast shocks because the tla/H^ 
ratios are too low. Extremely slow shocks (< 40 km s _1 ) 



M dMt (10 B M S ) 

Fig. 7. Mgas molecular hydrogen mass given by CO mea- 
surements as a function of M^ust obtained from IRAS 
data. The dust mass highly depends on the dust tempera- 
ture which is unknown : the dust temperature is estimated 
to 35K on the upper figure and to 40K on the figure be- 
low. The straight lines represent gas to dust ratio of 200 
(dashed), 500 (full) and 1500 (dotted). The temperature 
uncertainties exclude any clear conclusion about the cold 
gas origin. 



shocks, and are more consistent with the observations. 
But slow shocks are less efficient. The most likely solution 
is UV irradiation by very hot stars, implied by a star 
formation rate of only a few solar masses per year. A 
recent survey of H-band and K-band spectra in 32 central 
cluster galaxies have led to 23 detections in rovibrationnal 
H 2 lines, see Edge et al. (2002) and UV fluorescence 



excitation is ruled out. The molecular hydrogen is more 
probably thermally excited in dense gas with density 
exceeding 10 5 cm -3 at temperature evaluated between 
1000-2500 K. Young stars heating a population of dense 
clouds is invoked ( Wilma n et al. 2002|) . According to the 
authors, these dense regions might be self gravitating 
clouds deposited directly by the cooling flow or confined 
in high pressure behind strong shocks. Correlation 
are also shown between H a emission lines, warm H 2 
rovibrationnal lines and cold millimetric emission lines 
suggesting related exciting mechanisms of these different 
phases of the gas. The large masses of excited H 2 , around 
10 5_6 M Q could suggest that the cold molecular gas mass 
could have been underestimated (because of a lower 
metallicity for example) or is hidden in optically thick 
dense clouds, see Ferland et al. (2002) for a discussion 



of the physical conditions within dense cold clouds in 
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How much gas is deposited in cooling flows is still 
an open question. The gas cooling in the flow is probably 
multi-phase, and there are hints the CO detected here is 
the residual of the cooled gas. But this cold gas emission 
could also be due to subcluster structures, gas stripped 
from neighbouring galaxies or galactic clouds not seen 
until now and heated by mechanisms linked to the flow, 
like shocks or starburst. More investigations are required 
to explore the properties of this important component in 
cooling flow cluster cores. The study of the morphological 
structure of the cold gas and especially its dynamics will 
help to confirm its place in the flow. High resolution 
maps, obtained thanks to the 1RAM millimeter interfer- 
ometer, have been obtained for Abell 1795 in CO(l-O) 
and CO(2-l). These maps show an extended emission of 
the cold gas (Salome & Combes, in prep). They underline 
the possible link between the cold gas detected with the 
30m telescope and the cooling gas seen at higher energy. 
Recent OVRO observations by Edge & Frayer (2003) also 



show CO(l-O) emission maps in 5 cooling flow clusters 
of galaxies : A1068, RXJ0821+07, Zw3146, A1835 and 
RXJ0338+09. The authors conclude the gas previously 
detected with the single dish telescope is confined in 
the central region. More Plateau de Bure interferometric 
observations with higher sensitivity and spatial resolution 
are in progress now in RXJ0821+07 to see whether the 
cold gas is extended (as for Abell 1795) or centrally 
concentrated around the cD (as suggest the OVRO 
observations). Interferometric observations on a wider 
sample of CO detected cooling flow have now to be lead in 
order to explore the similarities and differences between 
clusters and definitively confirm the detection of the cold 
residual in cooling flows. 

5. Conclusions 

A sample of 32 cooling flow clusters of galaxies, selected 
on their mass deposition rate, and their Ha luminosity, 
have been observed in both CO(l-O) and CO(2-l) emis- 
sion lines. In total 6 clear detections are claimed, with 
4 other possible detections. Molecular hydrogen mass es- 
timates have been deduced for these galaxies and upper 
limits have been computed for the other ones. The derived 
MfHs) are up to 10 10 M Q in the 22 central arcseconds ob- 
served with the 30m telescope (that is typically the central 
^23kpc region at z=0.05). These masses appear to be re- 
lated to the cooling rate deduced from X-ray data : there 
is a trend of correlation with Mx results, and no longer 
large discrepancies between the mass deposition rates and 
the cold gas masses (according to recent mass deposition 
rates reevaluation from Chandra and XMM-Newton) . The 
apparent gas-to-dust ratio, derived from the CO emission 
and dust far-infrared emission is larger for the gas in cool- 
ing flow galaxies than in normal spirals, but uncertainties 
about the dust temperature preculdes any clear conclu- 
sions. The best correlation is between the cold gas masses 



(2001). Further work is to be done now to confirm that CO 
lines, revealed by single dish millimetric observations, are 
tracers of the long searched cold phase in cooling flows. In 
this context, more interferometric observations in CO(1-0) 
and CO(2-l) are required. 
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Fig. 8. CO(l-O) and CO(2-l) emission lines observed with the IRAM 30m telescope. On the Y-axis, main beam 
temperature (in mK) versus velocity (in km/s) on the X-axis . 
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